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Abstract 

One- and two-dimensional operators which originate from the asymptotic form of the three-body 
Coulomb wave equation in parabolic coordinates are treated within the context of square integrable 
basis set. The matrix representations of Green's functions corresponding to these operators are 
obtained. 
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I. INTRODUCTION 

The three Coulomb (3C) wave functions have been introduced to approximate the 
three-body continuum Coulomb wave function. The three-body Coulomb wave equation is 
asymptotically separable in terms of parabolic coordinates introduced by Klar 

= ris + ki, ■ Vis, rjj = ri, - kis ■ r^,, (1) 

where r^^ and k^^, /, s = 1, 2, 3 and / ^ s, are the relative coordinate and momentum vectors 
between the particles I and s. The 3C wave functions are proportional to solutions \E'3c of 
the asymptotic equation [3] (atomic units are used): 

where t/s = '^'f^^'° , uis = are the reduced masses. Here the operators hf, and h„. are 

defined by 

A solution \I^3(7 with pure outgoing behavior can be written as a product of three two-body 
Coulomb wave functions: 



with 



Uj{^j) = iFi{itis, 1; -ikis^j). (6) 

The solutions (JSj) are correct only when all particles are far from each other. In the last 
years many works have been published with different proposal to improve the simply 3C 
see and references therein). In this paper we seek to consider the possibilities 

for computing the three-body continuum Coulomb wave function which are afforded by 
expansion in a set of square-integrable functions. In principle, the short-range part of the 
Hamiltonian, which contains mixed second derivative terms in the expression for the kinetic 
energy [sl, can be approximated by a finite order matrix. The long-range part ([2]) of the 
Hamiltonian consists of two-dimensional operators 

i)j = h^^ + hr,^ + 2kistis + C^i^j +r]j), (7) 
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which are related by the constraint 

Ci + C2 + Cs = 0. (8) 

To treat the one-dimensional operators 

+ 2kt + C^, hr^ + 2kt + Ct] (9) 

within the context of square integrable basis set, the J-matrix method j3, or the tools of 
the "Tridiagonal Physics" program (see and reference therein) can be employed. In the 
framework of these methods a special basis set is used which supports an infinite tridiagonal 
matrix representation of the operator. Thus the resulting three-term recursion is analytically 
solved. 

In particular, in this paper we obtain matrices of Green's functions (resolvent) corre- 
sponding to the one-dimensional operators ([9]) and the two-dimensional operator 

\) = h^ + h^ + 2kto + C(e + v) (10) 

in the set of parabolic Sturmian functions, 

0nin2(^, V) = y^nliO ^n2{v), (H) 

ipnix) = V2be-'^L^{2bx), (12) 
where b is the scale parameter. This basis set have been used in the analysis of the Coulomb 



potential within the parabolic formulation of the J-matrix method 10| . The basis functions 
f lT^ are orthonormal: 

00 

dx Lfnix) Ifrnix) = Snm- (13) 





In Sec. II using the tridiagonal matrix representations of one-dimensional operators ([9]) in 
the bases (fT2l) . we construct the corresponding Green's function matrices. In this case we 
do not seek to determine the Green's matrices uniquely. In Sec. Ill the weight function 
is obtained for the orthogonal polynomials satisfying the three-term recurrence relation. 
The two-dimensional Green's function matrix elements are expressed as convolution of one- 
dimensional Green's matrix elements in Sec. III. An orthogonality relation employed in the 
two-dimensional Green's matrix construction is derived in the Appendix. 
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II. ONE-DIMENSIONAL COULOMB GREEN'S FUNCTION MATRICES 



a) C = 



The matrix representation + 2ktl^ [I^ is the unit matrix] of the operator + 2kt in 
the basis set {v^n(0}-^o P^j) is tridiagonal 

f bo d^ \ 
ai bi 3,2 

a2 &2 c?3 

03 X X 

^ XXX 



+ 2ktl(: 



(14) 



where 



(15) 



= (6 + ik) + 2bn + 2kt, 

o-n = {b — ik)n, dn = (b + ik)n. 
To construct the Green's function matrix g^, which is matrix inverse of the infinite tridiagonal 
matrix f|T^ . consider the three-term recurrence relation 

a„ Wn-l + bnWn + d^+l W„+l = 0, 71 > I. (16) 

It can be easily verified that 

(—^ \n v(r^ J- 1 — n-t\ 

(17) 



Vnit] = W Wl T-^ 2Fi{-n, It- -n + it- 0, 



V{l-it) 



where C, = is the "regular" solution of f|T6l) which satisfies the initial conditions: 

Po{t; = 1, bopoit; + dipiit; Q = 0. (18) 

Suffice it to say that apart from the factor (^2^)**5 Pn is the coefficient of the nth basis 
function fn{0 p2|) in the expansion of m(^) (EI). Note that pn are polynomials of degree n 
in t. 

The second solution g„ of the recursion ( |T6l) can be obtained from the condition that 
satisfies the same differential equation as the Pn [s]. In other words, ii Pn ~ 2Fi{a, b; c; z), 
then Qn ~ 2;^~'^2-Fi(cf — c + 1, 6 — c + 1; 2 — c; 2;). It is readily verified that an appropriate 
Quit; C) is 

9nit; C) = -^fgEf (-0"+^ 2^1(1 -tt,n + l;n + 2-tt; () 



n\r{l-it) I ( 
'r{n+2-it) lc-1 



n+1 



(19) 



iFi(n + l, n + 1; n + 2 - it; 



C-i' 



This satisfies the initial condition 



bo qoit; () + di qi{t; () = b - ik. 



(20) 



Multiplying [h^ + 2ktl^] by the diagonal matrix Z = [1, ( ^, . . . , ( " 
the symmetric tridiagonal matrix T: 



] , we obtain the 



T = Z [h^ + 2ktl^] 



"1 



tti Pi a2 

"2 ^2 "3 

as X X 

y XXX 



(21) 



with nonzero elements 



Pn bn/ C 5 '^n/ C 

Notice that p„ and g„ satisfy the three term recurrence relation 

anWn-l+ PnWn + an+lWn+l =0, U > 1. 



(22) 



(23) 



Thus, to invert the symmetric tridiagonal matrix T, one can draw on the standard method 
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12l |. Namely, the elements of a Green's matrix g-p, which is the matrix inverse to T, can 
be determined by 



(24) 



9nmW = w{(n^) ' ^ ^ mm(n, m), /i = max(n, m) 

where the Wronskian W is defined as 

W{q, p) = an [qn{t; C) Pn-i{t; C) - qn-i{t; OPn{t; 0] 

= b-ik = -2ik (^^) . 

The Green's matrix g^=[h^ + 2ktl^]~'^ is related to g^: g^ = grZ. Therefore, we can 
express the matrix g^ elements in the form 



(25) 



9nm (^) 



■Piy{t; Oq,i{t; C)- 



(26) 



2A; V C J C 

From ([3]) and (jlj) it follows that the Green's matrix g^=[h^ + 2fctl^]~^ and g^ are complex 
conjugates (for real k and t): 

a'ln.it) = {gUt))* . (27) 



Note that there is an ambiguity in determining the matrix g^, since the solution g„(t; C) 
is not unique: 

Ut; () = qn{t; C) + y{t)Pn{t] C), (28) 

where y{t) is an arbitrary function of t, also satisfies (fTB]l . 
b) C / 

It is not difficult to convince oneself that the differential equation 



is satisfied by the function 



u 



hi. + 2kt + Ci u{Ci = 



where 



X = ik 



1 _ 2C 1 



2t + i{l- Jl-^ 



7 = fc\/l-T^ 



2C 

fc2 



(29) 



(30) 



(31) 



Clearly, in the limit C = the function fl30l) reduces to IQ. 

Because the operator ^, evaluated in the basis (fT2|) . has the symmetric tri diagonal form 



'2b ^' 



Qn, n' ^ 



n' = n — 1, 
j-^{2n + l), n' = n, 
{n + 1), n' = n + 1, 



(32) 



the matrix representation of the operator + 2kt + is also tridiagonal. Thus, 

tti hi d2 

+ 2ktl^ + CQ^ = a2 62 4 

03 X X 

^ X X X y 



where 



(33) 



bn = {b+g + ik)+ 2(6 + g)n + 2kt, 
On = (^ - £ - ik)n, dn= {b- §■ + ik)n. 



(34) 
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Then, if we consider the coefficients of the regular solution (|30|) expansion in the basis 
set ( fT2l) . we obtain that 

snit; C)=e-pr,{T- 0, (35) 



where 



2C 



b - \ , c 

2b 

is the solution of the three term recurrence relation 



2C 



(36) 



a„ Wn-l + bn Wn + d^+l Wn+l =0, U > I. 

It is easy to verify that a second solution of Eq. flTTI) can be expressed in the form 

c„(t; C)=r+ig„(r; 0- 



(37) 



(38) 



Notice that the matrix + 2ktl^ + CQ^ inversion procedure is simplified if we introduce 
the symmetric tri diagonal matrix T: 

ai Pi oil 

012 ^2 as 
as X X 

\^ X X X y 



T = Z [h^ + 2kt\ + CQg] 



(39) 



Here Z is the diagonal matrix: Z=[l,x^,...,x'^,. 
and a„ are given by 

Pn bxij^X. 1 O^n dnlH. 

Then, it is easy to check that Sn and c„ also satisfy the recursion equation 



The elements /5„ 



a„ Wn-\ + /3„ U7„ + a„+i t(7„+i =0, > 1, 



(40) 



(41) 



and the Wronskian 



W{c,S)=an\Cn{t\ C)Sn-l{t- C)-C„_i(t; C)Sn{t- C)] 



(42) 



is independent of n. Namely, 



Finally, given the two linearly independent solutions s„ and c„ of Eq. (137|) . we can express 
the Green's matrix = [h^ + 2ktl^ + CQ^]^^ elements in the form 

C) g) ,-^ . 1. C). (44) 

Clearly, the Green's matrix g^=[h^ + 2A;tI^ + CQ^] ^ is obtained from by replacing 
t —t, k —k [this leaves kt unchanged]. 



III. ORTHOGONAL POLYNOMIALS p„ 

In this section we obtain the weight function with respect to which the polynomials p„ 
( ITTI) are orthonormal. 

The appropriate Kummer's relation [(15.3.7) in Ref. [l^] expresses the solution g„ as 

^ ^ (45) 



|arg(-C)| <vr. (46) 



Then introduce the function 



P«; C) ^ ^<^-'f'^'"-^'" (47) 

Zm 

and evaluate the integral 

J(C)= /" rftp(t; 0, (48) 



c 

where C runs along real axis, except for an infinitesimal indentation on the underside of the 
pole t = 0; see Fig. 1. 

To evaluate I{C) for \(\ > 1, change the variable of integration to s: t = transforming 
the path of integration into the circle Ci in the complex s-plane; see Fig. 2. Thus the interior 
of Ci corresponds to the upper half of the t-plane and hence the integral 

X(C) = / (jf^Pi(^5C), (49) 

Ci 

where p{t; () = p{^.] () = pi(s; (), is reduced to the sum of the residues at the poles 
[corresponding to t„ = in, n = 0, 1, . . .] inside the circle contour Ci. 



oo _ 

^(C)=^E7^' (50) 



n=0 



1. e. 



m = (51) 

In the case \C,\ < 1 it is convenient to change variables from t to s: t = The 
integration contour C is transformed into the circle C2 [see Fig. 3] and the lower half of the 
t-plane is mapped onto the interior of C2. Then introducing P2{s; C) = p(^; C) = pii\ C); 
we obtain 



{s + if^^''^' 1-C 

^ ' n=l 

Note that for /m(C) 7^ (gnD implies that 



C2 



/•It 

^ _ 27rf > < arg(C) < 0, 

p(t; = { (53) 



2^ — 7, < arg(C) < TT. 

Using the Sokhotsky formula we obtain 



m=vj dtpit- ()-'-. (54) 

—00 

In this case ( 153|) provides the convergence of the principal- value integral on the right hand 
side of (151. 

Clearly, integrals 

j dtpit; QPnit; QPnXt; C) (55) 
c 

can be expanded in terms of derivatives of (l48l) . It is now easy to obtain that 

^ (^^) / OPn(t; QPmit] C) = 5nm, (56) 

c 

i. e. p ( l47l) is the weight function for the orthogonal polynomials pn ( flTl) . 
For Im[C) 7^ we can rewrite the orthonormality relation (l56ll as 

^ (^) {^] ^^P^^'^ ^) " = (57) 

where the weight function p in the form fl53|) is used. Note that in the case C = we have 
C = e*''', < 0, and the weight function reduces to 

Po(^; C) = (58) 

It may be remarked that the polynomials p„ are discrete analogues of the charge parabolic 
Coulomb Sturmians introduced in [l^ . 



IV. TWO-DIMENSIONAL GREEN'S FUNCTION MATRICES 



In this section we obtain a matrix representation G{to; C) of the Green's function (re- 
solvent) of the operator flTUl) . Formally G(to; C) is the matrix inverse of the infinite matrix 

h = O + (g) + ® + C (Q^ ® + O Q^) , (59) 

i. e. 

hG(to; C) =1^®:^. (60) 

a) C7 = 

The elements of G(to) = G(to; 0) can be expressed as a convolution integral (see, e. g. 

Q) 

Cnma.mimal^o) = / dt g^^^^^t) g^^^^^to - t) , (61) 



Co 



where the integrand contains functions and g^^ that are proportional to the one- 
dimensional Green's function matrix elements [the non-uniqueness (1281) of the solution g„ is 
taken into account]: 

~9imit) = ^ (^) [A^qf^it; ()+x^{t)p^{t; 0]B^ (62) 

and 

9lJt) = ^iC - l)r Pui-t; C')[A,q,i-t; C~') + V.it) p,i-t; C')] D^. (63) 

Inserting fl6Tl) - fl63l) into fl60|) then gives the restriction on the factors and A^, functions 
{x„(t)}^Q and {yn{t)}^=Q, diagonal matrices B = [Bo,Bi, . . .] and D = [Do,Di, . . .], and 
the path of integration Co, namely, 

Co 

A(^B„i^ J dt gll^^^{tQ — t) = Vmin2^n2m2-, (65) 



Co 



Un-i 7712 ~l~ ^m-i 712 ■ (66) 



If we choose the contour C (see Fig. 1) as the path of integration in (16T|) [and (IMI) . 
(!65|l ]. we can draw on the orthonormality relation (!56|l to determine the rest of parameters. 
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In this case we readily check that, in particular, the set: = 0, Xn{t) = 2'Kipo{t; (), 
= Unit) ^ 0, Bn = -D„ = 1, Satisfies the conditions fl64l) - fl66l) . and hence the elements 
Gnin2,mim2 Can be expressed in the form 

Gmn2,mim2{to) = ^ (^) (vj dt Po{t; CKi(t; C) Prmit] C) (^0 " t) 

^ -~ . (67) 

Further, combining ( l56l) with ( lAll) and (1271) . we can obtain another allowable sets of the 
parameters. For instance, = 1, Xn{t) = 27iipo(t; (), Ar, = ^, yn{t) =0, -B„ = -D„ = 1. 



b) C / 

The elements of the matrix G(ito; 0) may be written as the convolution integral 

G'„in2, mi 012(^0; ^ J ^'^diimiit; C) gl!^^^^{to - t] C). (68) 



C3 



Here, the non- uniqueness of only is taken into consideration for simplicity, i. e. and 

are taken to be 

j / V — 1 \ 



^ A:2 / \ yl ^2 

The path of integration C3 in the complex r-plane (see Fig. 4) is identical to the contour C 
in the complex t-plane. The integration over r 0361) is performed on the assumption that r 
and C are independent of one another. 

Allowable parameters A^, x„(r) and Bm can be determined by substitution fl68l) - fl70|) 
into (!60|) . For instance, putting A^ = Q and x„(r) = 2A;p(r; (^), we obtain 3^. = 



Thus the matrix elements of the two-dimensional Green's function can 
be written in the form 



Grii 712, mi m2(t0i ^) 

/ X (71) 

= C^(¥) (^"'Jdrpir; OPniir; C)Pmi(r; C) ^7;;2 m2 (^o - t; C). 



C3 



11 



Appendix A: ONE USEFUL ORTHOGONALITY RELATION 



In this appendix we derive the orthogonahty relation 

dtgimif) = Sum 



2ik 

IT 



(Al) 



Using the integral representation [Eq. (15.3.1) in Ref. [13|] of the hypergeometric function 
in f|T9|) . we can rewrite the integral on the left hand side of flAll) in the form 

1 oo 

dx- — — I dtp,{t- ^)e-*"°(i--). (A2) 



c 



1 — x-^ 



vr vc-i; C" 

Notice that the integral over t in (lA2p consists of integrals 



[x-l) 



d 

dx 



5{x) 



(A3) 



i. e. is expressed in terms of derivatives of the delta function. Further, introducing the 
function 

VA^) = (A4) 



1 — x-^ 



we obtain 



Ije= / dx yj{x) S^^\x) 



l(-l)^yf(0) = f (-1)^+^ £ = /i, 



(A5) 



From ( 1A5I) it follows that the integral in (1A2I) is nonzero if n = m {= u = fi). In this case 
the only contribution to the integral comes from the leading term of the polynomial Pn{t', C) 
[which is equal to ^ (C — 1)" (—it)"]. Thus, we obtain for (1A2I) : 



(t^) ^^(C-1)" 27rX„„ = 1. 



If c 



(A6) 
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Figure 1: The path of integration C in (148p . 
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Figure 3: The integration contour C2 in (f52]l . 
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Figure 4: The path of integration C3 of the convolution integral in (|68i 
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